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ABSTRACT 

The terms of Sutcliffets development formula are evalu- 
ated by finite differences for one hundred cyclones. The 
values of relative divergence obtained are correlated with 
the subsequent 24-hour pressure change and 24-hour average 
speed. A correlation of -.08 was obtained between the relative 
divergence at the center and the subsequent 24-hour pressure 
change. Computing the algebraic sum of relative divergence 
across the center and correlating this value with the 24-hour 
pressure change gave a coefficient of .28. Using the alge- 
braic difference of relative divergence across the center 
and correlating this with the subsequent 24-hour average 
speed gave a coefficient of -.02. It therefore appears that 
the finite difference evaluation of the Sutcliffe's develop- 
ment equation, as made in this investigation, is not a useful 
method for a 24-hour prognosis of the deepening, filling, or 
movement of winter cyclones. 

The author wishes to express gratitude for the assistance 
of Professor G. J. Haltiner, Department of Aerology, U. 5S. 
Naval Postgraduate School, Monterey, California. 
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1. Introduction. 

Any attempt to forecast the future position or intensity 
of a cyclonic center should be based on the present dynamic 
processes as well as the changes which these processes will 
undergo. The latter part of this statement involves the more 
difficult phase of forecasting. There have been various 
approaches to the solution of this problem, none have been 
one hundred percent effective. In some instances the foree 
cast is entirely subjective and depends primarily on the 
experience of the forecaster. On the other hand, empirical 
or statistical rules have been devised which give satisfactory 
results in some geographical areas. It was the purpose of 
this investigation to see if Sutcliffe's development formula 

[5 | would serve as an accurate indication of the future in- 
tensity of a cyclonic center and the average speed it would 


maintain for a twenty-four hour period. 
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2. Resumé of the derivation of Sutcliffe's development 
formula. 


It may be desirable to briefly outline the derivation 
of this development formula and note the assumptions that 
are made. Essentially he begins with the general vorticity 


equation in the form 


Mf AV UplSeF) =~ CHW - bh (Ypex eh) 


st 


where pressure is taken as the vertical variable. According 
to atts nates by Sutcliffe, the last term, called the "tilt" 
term, is largely offset by 5 (Vp: VV) - Another modification 
is made by neglecting the vertical advection of vorticity 
which is nome l\y believed to be small. These simplifications 


reduce equation (1). to 


oh + V, * Vp (§+ F) = - F(Y, “V,) (2 ) 


Equation (2) involves the divergence at any one level, but 
the local field of pressure is not the result of the motion 
at one level, but rather of the integrated effects of the 
three dimensional motion. In order to get an estimate of 
the contribution of a layer, the difference is taken, result- 
ing in 

F(Vp°V, - Vp: Wyo) = 2 Get £) + Vue! Vp Gor f) -V, ° Vp (S+ F) @) 


where the subscript "0" denotes the lower layer. Applying 
the geostrophic approximation for vorticity in the form of 


-9 V? 
S=2 V,2 


2 





the local rate of change of the difference in vorticity 


becomes 


“2 3, (Vez -v2,) = - 2 YF dh y 


where "h" is the thickness between the two levels. Substi-e 


tution of equation (4) into equation (3) gives 


F(Vp-V,- We.) = - 2 V, 2b + Vy Up (8 +f) - Vu Ne (P+), (5) 


In order to use this equation'as a prognostic tool, the 
term of ab must be estimated or replaced by some reasonable 
assumption for the thickmess change. . In this derivation it 
{s assumed that temperature changes and therefore thickness 
changes are purely advective; and, finally, that this advection 
may be approximated by the geostrophic wind. That is 2b = -\V, Mh. 
With this assumption, noneadiabatic effects and vertical 
motions are neglected, and in some areas these may have marked 
affects on the thermal field. Using the advective assumption, 
as well as the equations, Vo =Vi + and 4 - Sgot f; 9 
the final form of equation (5) gives the "relative divergence" 


of Sutcliffe as 


‘Vio-T- ae. V 
(Ye Vi ~ Ve Vio) = we (2te + d+ dt) (le) 


(a) (b) (c) 
where s'! is the arc length in the direction of the thiclmess 


lines. 


We will now examine each of the various terms of this 


equatione 
&e From term (a) it is indicated that thermal wind over 
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a surface center of positive vorticity implies relative 
divergence ahead, relative convergence to the rear. However, 
this term has no effect on the intensification of the center, 
because at the center 23s is zero. Hence it merely affects 
the propagation of the wave. 

.-< aur This term indicates that there should be 
es ssiraen si where the thermal vorticity decreases 
in the direction of the thermal wind and relative convergence 
where the thermal vorticity increases in the direction of 
wind shear. In the normal synoptic situation, where the 
thickness trough is in the rear of a surface depression and 
a thicimess ridge over or ahead of it, there would be a de- 
crease of thermal vorticity downstream over the surface de- 
pression, which would mean relative divergence. On the other 
hand, when the thermal trough precedes the surface depression, 
there should be relative convergence. 

ce - af From this term one would get relative 
convergence when the thermal wind has a poleward component 
and relative divergence with an equatorward component. How- 
ever, in midelatitudes, the magnitude of the rate of change 
of the coriolis parameter is generally small compared to bse 
other terms and for purposes of this investigation will be 
neglected. 

It has been observed that divergence varies vertically 

and that the normal distribution over a surface low system 
is to have convergence in the lower layers with divergence 


above. In order to have development, the integrated divergence 
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aloft must over compensate for the integrated convergence in 
the lower layers. The amount of this over compensation 
should be a measure of the deepening or filling of the sur- 
face system. Instead of measuring the divergence and con- 
vergence for the entire vertical extent over a system, 
measuring it at some BpeciriSdasver or levels might be 
feasible and still give qualitatively similar results. 
Assuming a schematic vertical distribution of divergence 
over a cyclone or anticyclone as in Figure 1, Sutcliffe 
implies that the differences of divergence of any two levels 
should be an indication of the deepening or filling of the 


surface system. 





cony, ° diy Con V. ody. 
Simple Cyclonic Type Simple Anticyconic Type 
Figure 1 


Schematic representation of the distribution of divergence. 


'That is, if the relative divergence over a cyclone is a large 
positive value, the deepening should be large and as this 
value decreases, the deepening should also decrease. An 
equivalent statement for anticyclones would associate negative 


values of relative divergence with anticyclonic development. 
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5. Investigation procedures and correlations. 

The purpose of this investigation was to see mw it 
would be possible to correlate the present values of relative 
divergence as measured by terms (a) and (b) of equation (6) 
with a twenty-four hour deepening or filling of a surface 
low, as well as the average speed of Penonsnt over the same 
period of time. One hundred cases of surface lows in the 
period of December 1956 through February 1957 were used for 
the investigation. These lows were located in mid-latitudes, 
ranging from 28°N to 60°N. They were selected primarily over 
the North American continent in order to have the advantage 
of a dense reporting network, thereby reducing the analysis 
error due to subjectivity. 

The initial attempt in the investigation was to evaluate 
the relative divergence at the center of a surface low. In 
regard to Sutoliffe's formula this meant evaluating the term 
aa ar » The other term would have no contribution because 
the surface vorticity center and low center are normally co- 
incident, thereby making 2 fo equal to zero. Positions 
three hundred nautical miles upstream and downstream along 
the same thickmess line were used and the calculation of 


thermal vorticity at these points was made by the finite 


adifference approximation of 2 Vih » namely, 
1 (a4) «- A grid distance, dad, of 5 degrees latitude 


was used which is consistent with the values used in graphical 
and numerical techniques. A second finite difference approx- 


imation was used to evaluate the rate of change of thermal 
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vorticity across the center. The latitude of the surface 
low center and the geostrophic thermal wind over the center 
were the other values used in the evaluation. 

A plot of the relative divergence, arrived at by this 
method, versus the twenty-four hour pressure change is shown 
in Figure 2. There is no definite grouping of the values, 
which indicates that the future deepening or filling of the 
system cannot be predicted using the relative divergence at 
the center. A correlation coefficient of -.08 was computed 
from the data of this evaluation. 

The results of the initial investigation were not cone 
clusive so another approach was taken. In accordance with 
the kinematic formula of deepening or filling of the surface 
center should be proportional to the algebraic sum of the 
divergence ahead and behind the system C1) « That Legman 
there is more relative divergence ahead of the system than 
there is relative convergence behind the system, the system 
should deepen. To test this, it was necessary to evaluate 
terms (a) and (b) of Sutcliffe'ts development formula. The 
values of thermal vorticity and surface vorticity for positions 
along the same thickness line six hundred nautical miles ahead 
and behind the surface low center, as well as at the center 
were calculated. The geostrophic finite difference approxi- 
mations of z Vir and 4 Vz. were used. 
The finite difference approximation was used to compute the 
rate of change of these quantities along the thickness lines 
for points three hundred nautical miles ahead and behind the 
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center. The thermal wind and latitude of the points three 
hundred nautical miles ahead and behind the surface center 
were the other values used in the computation. The algebraic 
sum of terms (a) and (b) were computed across the low and 
plotted versus the actual deepening or filling for twenty- 
four hours. Figure 3 shows this plot and the scatter of 
this diagram implies that very little can be inferred from 
the present values of relative divergence sums and the subse- 
quent pressure change. A correlation coefficient of .28 
between the relative divergence sums and the twenty-four 
hour pressure change was computed from the data in this test. 
The third phase of the investigation consisted of an 
attempt to correlate the relative divergence with the twenty- 
four hour speed of the surface low system. The speed of the 
system should be proportional to the gradient of relative 
divergence across the system. To evaluate this quantity, 
the algebraic difference of the relative divergence values 
for the points three hundred nautical miles ahead and behind 
was taken. The subsequent speed of the associated low was 
evaluated in units of degrees of latitude per twenty-four 
hourse A scatter diagram of the speed versus the difference 
of relative divergence is shown in Figure 4. The correlation 
coefficient of these two quantities is -.02, and thus no 
definite statement can be made as to what average speed can 
be expected from the value of the difference of relative 


divergence across the system. 
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4. Results and interpretations. 

As a result of the observations and computations made 
in this investigation, it appears that on the basis of 
Sutcliffe's development formula the future development or 
displacement cannot be forecast. 

The poor relationships resulting from the various phases 
of the investigation soula be attributed to the fact that 
the actual vertical distribution of divergence, as well as 
relative divergence is more complex than that which was 
assumed by Sutcliffe. Also some of the assumptions ine 
volved in the derivation of the development formula are not 
completely realized at all times. Perhaps the length of 
time that was chosen for the future intensity was too long, 
but on a cursory examination of the twelve-hour pressure 
changes and speed, no improvement was apparent. It also 
may be indicated that the distances used in the evaluation 
may be too large and that positions closer to the center 
would give better approximations for the finite difference 
method. Perhaps making some modifications in the method of 
evaluation, or by empirically allowing for the neglected 
terms of the development formula, a more promising result 
could be obtained. An evaluation of the terms of Sutcliffe's 
development formula on a synoptic chart by use of a small 
grid system and Taylor Series expansion has been carried out 
by Sawyer and Matthewman (el, but no attempt was indicated 
of correlating the results with the future intensity of the 


~ 


cyclonic system. 
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